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Monogenic functions in the algebra of 5-dimensional spacetime have been used previously by the 
author as first principle in different areas of fundamental physics; the paper recovers that principle 
applying it to the hydrogen atom. The equation that results from the monogenic condition is 
formally equivalent to Dirac's and so its solutions resemble closely those found in the literature. 
The use of the monogenic condition as point of departure as not only the advantage of being a 
unified approach but also provides very strong links with geometry that are completely lost in the 
usual approach. 



PACS numbers: 12.20.-m; ll.lO.Kk. 



I. INTRODUCTION 

I have been advocating in recent papers that the ma- 
jority of physics equations can be derived from an ap- 
propriately chosen geometry by exploration of the mono- 
genic condition. Monogenic functions are not familiar to 
everybody but they are really the natural extension of 
analytic functions when one uses the formalism of geo- 
metric algebra [U, [3, S| ; those functions zero the vector 
derivative defined on the geometric algebra of the partic- 
ular geometry under study. 

In I showed how special relativity and the Dirac 
equation could be derived from the monogenic condition 
applied in the geometric algebra of 5-dimensional space- 
time G4_i. An earlier paper Q proved that the same 
condition in the same algebra was sufficient to produce a 
symmetry group isomorphic to the standard model gauge 
group; unfortunately this paper is incorrect in the for- 
mulation of particle dynamics but the flaw was recently 
corrected; @ the same work introduces electrodynamics 
and electromagnetism in the monogenic formalism. Cos- 
mological consequences were drawn from the addition of 
an hyperspherical symmetry hypothesis with the conse- 
quent choice of hyperspherical coordinates. Summing 
up all those cited papers, I wrote a long book chapter, [g 

The present paper uses the 5D monogenic condition to 
study electron orbitals around a positively charged point 
sized nucleus, a problem known in quantum mechanics as 
the "Hydrogen atom problem." The equation that arises 
from the monogenic condition applied in 6*4,1 algebra is 
formally equivalent to the Dirac equation and so one ex- 
pects that its solutions are also formally equivalent to 
those obtained in relativistic quantum mechanics. It will 
be shown, however, that obtaining the solutions is greatly 
facilitated by the use of geometric algebra formalism and 
also that the resulting formulas are much more compact 
than standard ones and lend themselves to an easier ge- 
ometrical interpretation. The derivations make use of 
some methods and strategies drawn from Refs. P, Q but 



depart from those works in many important aspects. 



II. SOME GEOMETRIC ALGEBRA 

Geometric algebra is not usually taught in university 
courses and its presence in the literature is scarce; good 
reference works are (T], [3, S|- We will concentrate on 
the algebra of 5-dimensional spacetime because this will 
be our main working space; this algebra incorporates as 
subalgebras those of the usual 3-dimensional Euclidean 
space, Euclidean 4-space and Minkowski spacetime. We 
begin with the simpler 5-D flat space and progress to a 
5-D spacetime of general curvature (see Appendix [C] for 
more details.) 

The geometric algebra 6*4,1 of the hyperbolic 5- 
dimensional space we consider is generated by the co- 
ordinate frame of orthonormal basis vectors CTq such that 



iaof = -1, 
i'^if = 1, 

i^a'f^fj — 0, a ^ (3. 



(1) 
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Note that the English characters i, j, k range from 1 to 4 
while the Greek characters a,/?, 7 range from to 4. See 
Appendix |X] for the complete notation convention used. 

Any two basis vectors can be multiplied, producing 
the new entity called a bivector. This bivector is the ge- 
ometric product or, quite simply, the product, and it is 
distributive. Similarly to the product of two basis vec- 
tors, the product of three different basis vectors produces 
a trivector and so forth up to the fivevector, because five 
is the dimension of space. 

We will simplify the notation for basis vector products 
using multiple indices, i.e. UaOp = Uap- The algebra is 
32-dimensional and is spanned by the basis 

• 1 scalar, 1, 

• 5 vectors, (Tq, 

• 10 bivectors (area), ctq^, 

• 10 trivectors (volume), CTq/j^, 
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• 5 tetravectors (4- volume), io-Q, 

• 1 pseudoscalar (5- volume), i = (Toi234- 
Several elements of this basis square to unity: 

The remaining basis elements square to —1: 



(2) 



(^o)' = K)2 - {a,,k? = (i^.)' = = -1- (3) 

Note that the pseudoscalar i commutes with all the other 
basis elements while being a square root of — 1; this makes 
it a very special element which can play the role of the 
scalar imaginary in complex algebra. 

In 5-dimensional spacetime of general curvature, 
spanned by 5 coordinate frame vectors ga, the indices 
follow the conventions set forth in Appendix |^ We will 
also assume this spacetime to be a metric space whose 
metric tensor is given by 



(4) 



the double index is used with g to denote the inner prod- 
uct of frame vectors and not their geometric product. 
The space signature is ( — h + + +), which amounts to 
saying that t/oo < and ga > 0. If the coordinate frame 
vectors can be expressed as a linear combination of the 
orthonormed ones, we have 
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(5) 



where n^a is called the refractive index tensor or simply 
the refractive index; its 25 elements can vary from point 
to point as a function of the coordinates. 3 In this work 
we will not consider spaces of general curvature but only 
those verifying condition ([5]); in those spaces we define 
the vector and covariant derivatives (see appendix [Pjl. 



III. DIRAC'S EQUATION 

There is a class of functions of great importance, called 
monogenic functions,^ characterized by having null vec- 
tor derivative; a function is monogenic in flat space if 
and only if 



V^- = 0. 



(6) 



A monogenic function is not usually a scalar and has by 
necessity null Laplacian, as can be seen by dotting Eq. 
(|n]) with V on the left; it has solutions of the type 



(7) 



The Dirac equation can be derived from the monogenic 
condition, as shown in Q and briefly remembered here. 
For this effect it will be convenient to expand the mono- 
genic condition © as ( +a'^d4)'if = 0. If this is applied 



to the solution ([7]) and the derivative with respect to x'^ 
is evaluated we get 



(8) 



Let us now multiply both sides of the equation on the left 
by and note that bivector cr"*" squares to the identity 
while the 3 bivectors cr'*'" square to minus identity; we 
rename these bivectors as 7-bivectors in the form 7^ = 



a f^. Rewriting the equation in this form we get 
(7^a^ + ip4)*-0. 



(9) 



Invoking the isomorphism between G^^i and the complex 
algebra of 4 * 4 matrices, the only thing this equation 
needs to be recognized as Dirac's is the replacement of 
Pi by the particle's mass m; simultaneously we assign 
the energy E to po and 3D momentum p to cr'"pm . [J, [y] 
Alternatively we can multiply both sides of Eq. ([5]) by a" 
on the left to obtain the a, f3 form of the Dirac equation, 
through the assignment a™ = cr™°, /? = cr^°. 

Applying the monogenic condition to Eq. ([7]) we see 
that the following equation must be verified 

*o(a>„) = 0. (10) 

It has been shown Q that p = a^pa is a null vector and 
the wavefunction in Eq. ([7]) can be given a different form, 
taking in consideration the previous assignments 



= A{±aim + p + cro-E)e 



,i(Bt+p-x±mT+(5). 



(11) 



where A is the amplitude and x the 3- 

dimensional position. Although 'J does not look like a 
column Dirac spinor, it has the same number of compo- 
nents and can be written in that form, if desired. 

In Eq. ([8]) we used i as imaginary in the exponent but 
the monogenic condition would be equally verified if we 
had chosen any other algebra element whose square was 
minus unity; in Q the various possible such elements are 
analysed and discussed; in the present work we will not 
explore such possibility. Applying the vector derivative 
to this solution we have 

V* = (£;ct'' ± TOcr'')i* + V* = 0. (12) 



IV. ELECTROMAGNETIC POTENTIALS 

When working in curved spaces the monogenic condi- 
tion is naturally modified, replacing the vector derivative 
V with the covariant derivative D. A generalized mono- 
genic function is then a function that verifies the equation 



D^f = 0. 



(13) 



Remembering that an electron has minus unit charge 
in our units' system, for an electron with rest mass m in 
the presence of an electromagnetic potential A we have 
to consider the reciprocal frame (see @ and [3]) 



-A., 



(14) 
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which corresponds to the refractive index frame 

90 = o-f, + — 0-4, 94 = 0-4; (15) 
m 

it is easily verified that 9"'9f3 — 5"' p, as required by the 
definition of reciprocal frame. The covariant derivative 
is then 

D = cT°ao + V+ f-^a'^ + d'*^ 94. (16) 
\ m J 

We expect solutions of Eq. that are harmonic in t 
and T, which we write as 4* = '^(x'") exp[i(_Bt — ttit)], 
selecting the signs in the exponent for a forward propa- 
gating wave. The monogenic condition becomes 

[a^{E - Ao) - Vi - cr"M„, - a^m] * = 0. (17) 

For the following derivations we will follow closely the 
procedure explained in detail in Doran and Lasenby 
making the necessary changes to conform to our mono- 
genic formalism. In the case of a central field we have to 
make Am = and Aq = 4>(r), with r the radial coordi- 
nate and the previous equation becomes 

[(T°(£;- 0) - Vi- crV] * = 0. (18) 

Under those conditions wc write the previous equation 
as 

V* = P*, 

V -{Ea° - (jjcj" - ma% (19) 

The wavefunction ^' is multiplied by a vector on both 
sides of Eq. (fTQ]) but the vector on the Ihs has only 3D 
components while the vector on the rhs has only a'^ and 
tJ^ components; this arrangement is particularly useful to 
study commutativity of operators, as we shall see. Sup- 
pose we have an operator O that commutes with both V 
and 2? and suppose also that "if is in an eigenstate of O, 
that is = o^f; Eq. is then automatically verified. 
Our task is then to find the operators that verify the 
commutation conditions and solve for their eigenstates; 
this procedure is very similar to what is done in quantum 
mechanics for the hydrogen atom. 

In Appendix |E] we define the total angular momen- 
tum J3 whose eigenstates define the angular solutions; 
another associated operator is 

IC ^ (1 - xAV)a^° . (20) 

/C anti-commutes with V because the bivectors present 
in xAV belong to 3D space and a'^° anti-commutes with 
both cr" and ti^. The fact that it does not commute 
with V is not a problem as long as /C has symmetric 
eigenvalues; it commutes with V as shown in Appendix 
IeI So, we can assume to be in an eigenstate of K,, that 
is 

/C* = fc*. (21) 



With central potentials we will use spherical coordi- 
nates; this implies that the frame vectors undergo rota- 
tions and their derivatives must be considered as follows 

deOB = -cTr, d^ao = cos^ct^, (22) 
dga^ = 0, d^a^p — — sin 9 Or — cos 9 ag . 

The 3D part of the vector derivative becomes 

V = (7"^ dr + - {cr'^dg + CSC 9 cr'^d^). (23) 



V. ANGULAR SOLUTIONS 

In order to solve Eq. we make the ansatz "if = 
RYf , where y/ contains all the angular dependence and 
i? is a function of r, t and r only. We note also that 
xA Vi? = 0; in terms of this operator we can write 

xAV* = (1 -CT^'^fc)* = i?xAVy/. (24) 

The analysis starts with 3D monogenic functions, or 
spherical harmonics, defined by Vi/'t = 0; these functions 
are of the type = r*Yf{d, ip) and the vector derivative 
is 

VV't = tcrV*-ir/ + r* Vy/. (25) 
Since the first member must be null 

Vy/ = --a-y/. (26) 
r 

In terms of operator xA V the equation is 

-xAVF/ = ir/; (27) 

comparing with Eq. ([M]) we see that k = t + 1. Equa- 
tion (|27p is an eigenvalue equation satisfied by spherical 
harmonics with the general formula [H 

= [(s + t + i)p/(cos0) 

-P/+i(cos6l)a''^] e^'^'^", (28) 

where Pi{x) are Legendre polynomials, i > and — 1 — 
t < s < t; solutions for t < ~2 can be found through the 
relation 

-xAV((t''V'(t3) = -(t + 2)CT''V'CT^ (29) 

Because k = t + 1, the allowed values for k are given 
by |fc| > 0; this ensures that the eigenvalues of K, opera- 
tor are indeed symmetric, as we required in the previous 
section. 

A particularly simple formula for 3D monogenic func- 
tions is obtained when the two quantum numbers are 
identical, in which case we obtain 

V't = (rsin6')*e-*'^'"''. (30) 

Since Eq. (pS)) does not depend on s, we can always use 
this particular case for the discussion of radial solutions. 
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VI. RADIAL SOLUTIONS 

For an hydrogen-like atom we make (jj = Za/r, with Z 
the positive charge of the nucleus and a the fine structure 
constant. We will multiply Eq. on the left by x 
noting that xV — rdr + xA V; considering Eq. (|24p 



Za 



rdr-^ + (1 - fccr*")* = r i^Ea'''' - ^ a'"^ - ma'''"j i^-. 

(31) 

Introducing the function T — rVl/, this can be rearranged 
isolating drT 

drT = (Ecr°'- + md^niT + - (ka*° + Zaa''^\)T. (32) 

r 

It is useful to define the two multivectors 



F = -(£;o-"'' -mcr-*'')i, 
G = -{ka^^ + Zaa'-^i)- 

so that the monogenic condition becomes 



9,.T + F 



G 



The F and G operators satisfy 

G2 = fc2 - {Zaf = 



(33) 



(34) 



(35) 



They also observe the anticommutation relation 

FG + GF = ~2ZaE. (36) 

We will now make a change of variable to allow the 
separation of large and small r behaviour; this is y = fr 
with which we write 

T = $e-^. (37) 
The new function <I> satisfies 

+ ^ $ + (^y - 1^ * = 0. (38) 

We can certainly express $ as a power series; moreover 
this power series must not be infinite, otherwise it would 
not fall to zero at large r. Calling C; to the series coeffi- 
cients 



(39) 



The coefficients verify the recursion relation 



F 



{l + s + G)Ci = - [j -l]Ci.i. (40) 



For ; = it is 



Multiplying on the left by (s — G) we can see that = 
G^ = ly^ and we set s = v to avoid a central singularity. 

Since the series terminates at Z = n, the coefficient 
Cn+i must be null but it must still verify the recursion 
relation 



1 C„ = 0, 



(42) 



and so FC„ = /C„. Multiplying both sides of the recur- 
sion relation on the left by (F/f + 1) and replacing s by 
v 



f 



+ l]{n + i^ + G)Gn = 



F 



1 



F 



J J \.f 
Combining the two equations we get 



1 C„- 



0, 



2{n + iy) + G+jG 



This can in turn be manipulated to give 



and finally 



2{n + „) + i{GF + EG) 



ZaE 
n + V — = 0. 



0; 



/ 



(43) 



(44) 



(45) 



(46) 



This is the energy quantization equation, which can be 
arranged into the usual form by first manipulating to 



E 



n + V 



^/m? - E'^ Za 
and then rearranging to 

(Zaf 



E^ = 



(47) 



(48) 



VII. GROUND STATE OF THE HYDROGEN 
ATOM 

It is useful to analyse the ground state in order to 
find the character of the series coefficients. Since this 
is a spherically symmetric solution we ignore the 9 and 
Lp dependence. At large r we can neglect the angular 
dependence of cr'" and a second order equation can be 
written for iIj{x™') as 

drr^ « {Ea" - ma^) a'^dri^i = (m^ - E'^)^^. (49) 

Since we are looking for bound states it must he E < m 
and ij} goes with exp(— /r), with / = V m? — E'^ ; the 
large r solution is 



(s + G)Go = 0. 



(41) 



■0 = ipie . 



(50) 
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Inserting into Eq. ((T3|) 



(51) 



The equation is solved if ■i/'i contains a factor {Ea'^ 



For small r Eq. becomes 

V'0 « a^ipi. 

r 



(52) 



Because we are assuming -0 to be a radial function we try 
a solution of the type 



ij} = aa'' + bcr", 



(53) 



where a and b are scalar functions of r. We note that 
Va'' — 2/r and insert in Eq. ((52)l to get 



a' + — + b'a'-° = -^a^aa'' + 6a")i 
r ' 

= (a(7'-o + fe)fii. 
This implies the simultaneous equations 
, 2a Zab , 



Zaa . 
b — 1 



(54) 



(55) 



From the first equation in this set we take the derivative 
of b as 



b' = -{3a' + ra )— , 
Za 



(56) 



We can now combine the two equations for b' into the 
single differential equation 



Zaa 



3a' + ra 
Za 



which we solve as 



with 



V 



-1±^1 -{Zay 



(57) 

(58) 
(59) 



we will later argue that only the plus sign is physically 
meaningful. Inserting a into the second Eq. (|55p and 
solving we get for b 



b^^r\ 



(60) 



Actually it is convenient to multiply a and b by i. Sum- 
marizing the results for small and large r, respectively, 



Za 



■a + a 1 r 



r small, 



I V 1 



cr'i^ e r large. 



(61) 



In order to make the two factors in brackets compatible 
we must include a factor (1 — cr^") to get 



(62) 



I V / 



Za 



a' i (1 - a^")e 



40^„-/r• 



— (aO - a^) + (1 - a«)a'^i 



E + m 
~1~ 



(a" - a^) + (1 - a^)^^ 



(63) 



-fr 



The complete radial solution for the ground state is 
Za 



-77 



(a° - a^) + (1 - a«)a'-i 



r''e"■'^^ 



under the condition 

E + m Za 

./ " ^ 1 T v^l - (Za)2 

The left hand side can be written as 

E + m f I'm j 1'^ 



f 1 - E/m 1 ^ ^1 - (//m)2 



(64) 



(65) 



(66) 



In order to have a positive energy it must be / = mZa 
and 



E = m^/l - (Za)2 



(67) 



This is a particular case of the general energy quantiza- 
tion formula derived above, for n = and k = 1. 



VIII. CONCLUSION 

The relativistic solutions for the hydrogen atom can be 
derived from the monogenic condition applied to func- 
tions in the algebra of 5D spacetime. When the space is 
bent by the consideration of a central potential, the equa- 
tion that follow from the monogenic condition is formally 
equivalent to Dirac's and so one could expect equivalent 
solutions from the onset. The paper revisits the solutions 
of that equation applying the formalism consistent with 
the monogenic condition. 

There is nothing fundamentally new in the paper, in 
the sense that the energy levels that one obtains with 
the monogenic formalism are the same that can be found 
in the literature. However, the monogenic condition had 
previously been applied in other areas of physics showing 
high unifying potential; this work is just one further step 
in the path of unification. Besides that, the author be- 
lieves the monogenic formalism to be easier to apprehend 
than the more usual matrix formalism; this is a question 
of taste, though. 
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APPENDIX A: INDEXING CONVENTIONS 

In this section we establish the indexing conventions 
used in the paper. We deal with 5-dimensional space 
but we are also interested in two of its 4-dimensional 
subspaces and one 3-dimensional subspace; ideally our 
choice of indices should clearly identify their ranges in 
order to avoid the need to specify the latter in every 
equation. The diagram in Fig. [1] shows the index naming 
convention used in this paper; Einstein's summation con- 



a, 7 - 



{01,2;}, f} 



FIG. 1: Indices in the range {0, 4} will be denoted with Greek 
letters q, /3, 7. Indices in the range {0, 3} will also receive 
Greek letters but chosen from \i,v,(^. For indices in the range 
{1,4} we will use Latin letters i,j,k and finally for indices 
in the range {1,3} we will use also Latin letters chosen from 
m, n, o. 

vention will be adopted as well as the compact notation 
for partial derivatives da — d/dx"". 



APPENDIX B: NON-DIMENSIONAL UNITS 

The interpretation of t and t as time coordinates im- 
plies the use of a scale parameter which is naturally cho- 
sen as the vacuum speed of light c. We don't need to 
include this constant in our equations because we can 
always recover time intervals, if needed, introducing the 
speed of light at a later stage. We can even go a step fur- 
ther and eliminate all units from our equations so that 
they become pure number equations; in this way we will 
avoid cumbersome constants whenever coordinates have 
to appear as arguments of exponentials or trigonomet- 
ric functions. We note that, at least for the macroscopic 
world, physical units can all be reduced to four funda- 
mental ones; we can, for instance, choose length, time, 
mass and electric charge as fundamental, as we could 
just as well have chosen others. Measurements are then 
made by comparison with standards; of course we need 
four standards, one for each fundamental unit. But now 
note that there are four fundamental constants: Planck 
constant {H), gravitational constant (G), speed of light 
in vacuum (c) and proton electric charge (e), with which 
we can build four standards for the fundamental units. 
Table [T] lists the standards of this units' system, fre- 
quently called Planck units, which the authors prefer to 
designate by non-dimensional units. In this system all 
the fundamental constants, h, G, c, e, become unity, a 
particle's Compton frequency, defined hy v = mc^ /h, be- 
comes equal to the particle's mass and the frequent term 
GM/{c^r) is simplified to M/r. We can, in fact, take 



TABLE L Standards for non-dimensional units' system 
Length Time Mass Charge 



all measures to be non-dimensional, since the standards 
are defined with recourse to universal constants; this will 
be our posture. Geometry and physics become relations 
between pure numbers, vectors, bivectors, etc. and the 
geometric concept of distance is needed only for graphical 
representation. 



APPENDIX C: SOME COMPLEMENTS OF 
GEOMETRIC ALGEBRA 

In this section we expand the concepts given in Sec.lTTl 
introducing some useful relations and definitions. Start- 
ing with the basis elements that square to unity Eq. 
repeated here. 



(CI) 



it is easy to verify any of the above equations; sup- 
pose we want to check that (croy)^ = 1- Start by ex- 
panding the square and remove the compact notation 
(foy )^ = o'o'^i'^jO'oO'iaj, then swap the last aj twice to 
bring it next to its homonymous; each swap changes the 
sign, so an even number of swaps preserves the sign: 
(foy)^ — <^oai{aj)^aoai. From the third equation ([T|) 
we know that the squared vector is unity and we get suc- 
cessively ((Toy)^ = (JQaiaoo-i = -((To)^(cri)^ = -(o-q)^; 
using the first equation ([T]) we get finally (troy )^ = 1 as 
desired. 

The remaining basis elements square to —1 as can be 
verified in a similar manner, Eq. ([3]): 



(C2) 



Note that the pseudoscalar i commutes with all the other 
basis elements while being a square root of —1; this makes 
it a very special element which can play the role of the 
scalar imaginary in complex algebra. 

We can now address the geometric product of any two 
vectors a = a'^Ua and h = b^ap making use of the dis- 
tributive property 



ab - 



I -a"6" + ^'b' j -t- a"b<^aap; (C3) 



and we notice it can be decomposed into a symmetric 
part, a scalar called the inner or interior product, and 
an anti-symmetric part, a bivector called the outer or 
exterior product. 



ab = a-b + aAb, ba = a-b — aAb. 



(C4) 
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Reversing the definition one can write inner and outer 
products as 



a-b= - {ah + ba), 



a^b=^{ah-ha). (C5) 



The inner product is the same as the usual "dot prod- 
uct," the only difference being in the negative sign of 
the 0060 term; this is to be expected and is similar to 
what one finds in special relativity. The outer prod- 
uct represents an oriented area; in Euclidean 3-spacc 
it can be linked to the "cross product" by the relation 
cross(a, b) = — cri23aAb; here we introduced bold char- 
acters for 3-dimensional vectors and avoided defining a 
symbol for the cross product because we will not use it 
again. We also used the convention that interior and ex- 
terior products take precedence over geometric product 
in an expression. 

When a vector is operated with a multivector the inner 
product reduces the grade of each element by one unit 
and the outer product increases the grade by one. We 
will generalize the definition of inner and outer products 
below; under this generalized definition the inner product 
between a vector and a scalar produces a vector. Given 
a multivector a we refer to its gradc-r part by writing 
<a>r] the scalar or grade zero part is simply designated 
as < a >. By operating a vector with itself we obtain a 
scalar equal to the square of the vector's length 



aa = a-a + aAa = a-a. 



(C6) 



The definitions of inner and outer products can be ex- 
tended to general multivectors 



a-b = ^ (<a>a <b>i3) 



\a-D\ 



(C7) 



a,/3 



aA6 = ^(<a>a <6>/3>„+^. (C8) 

Two other useful products are the scalar product, de- 
noted as <ab> and commutator product, defined by 



a X 6 = - {ah — ha). 



(C9) 



In mixed product expressions we will use the convention 
that inner and outer products take precedence over geo- 
metric products as said above. 

We will encounter exponentials with multivector expo- 
nents; two particular cases of exponentiation are specially 
important. If u is such that v? = —1 and ^ is a scalar 



Conversely if h is such that h =1 



jie 



= i + he + — + h 



3! 4! 



+ ... 



+h6 + h— + ...{=h sinh 6} 
= cosh 6 + h sinh 6. 



(Cll) 



The exponential of bivcctors is useful for defining rota- 
tions; a rotation of vector a by angle 6 on the ai2 plane 
is performed by 

a' = e^^'^/^ae"''^/^ = RaR; (C12) 

the tilde denotes reversion and reverses the order of all 
products. As a check we make a = ai 



cos ^ - (712 sin - ) ai 



* \ 2 ~^ 2/ ^^^"^^ 

= cos 0ai + sin 9a2 ■ 

Similarly, if we had made a = a2, the result would have 
been — sin^ai -|- cos 0a2- 

If we use B to represent a bivector whose plane is nor- 
mal to fJo and define its norm by|S| = a general 
rotation in 4-space is represented by the rotor 



R = 



cos 



\B\ 



B 



TT^ sm , 
LB V 2 



\B\ 



(CM) 



The rotation angle is |_B| and the rotation plane is defined 
by B. A rotor is defined as a unitary even multivector 
(a multivector with even grade components only) which 
squares to unity; we are particularly interested in rotors 
with bivector components. It is more general to define 
a rotation by a plane (bivector) then by an axis (vector) 
because the latter only works in 3D while the former is 
applicable in any dimension. When the plane of bivector 
B contains cto, a similar operation docs not produce a 
rotation but produces a boost instead. Take for instance 
B = (701^/2 and define the transformation operator T = 
exp(B); a transformation of the basis vector ao produces 
a' = f anT = e-'^^i^/^ e'^"^^/^ 



cosh ^ - (701 sinh | ) ctq 



* ^cosh — + (7qi sinh — 
= cosh 9ao + sinh 9ai . 



(C15) 



APPENDIX D: RECIPROCAL FRAME AND 
DERIVATIVE OPERATORS 



+u9~u— + ...{=usm9} (CIO) 



= cosf-t-wsmf. 



A reciprocal frame is defined by the condition 

5"-5/3 = <5"/3. (Dl) 
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Definine; g"P as the inverse of ga[3^ the matrix product of 
the two must be the identity matrix, which we can state 
as 

g'^-'gp^ = 5'^^. (D2) 
Using the definition ^ we have 

{9^^9^)-9p^5%- (D3) 
comparing with Eq. (|Dip we determine g"' 

It would be easy to verify that it is also g"^^ = 9" ■ 9^ 
and ga = 5075'^- 

In many situations of great interest the frame vectors 
ga can be expressed in terms of an orthonormed frame 
given by Eqs. (P). If the frame vectors can be expressed 
as linear combination of the orthonormed ones we have 

9a=n'^acrp, (D5) 

where n^a is called the refractive index tensor or simply 
the refractive index as said in the main text. When the 
refractive index is the identity we have ga = <Jol for the 
main or direct frame and g^ = —uq, g^ = Oi for the 
reciprocal frame, so that Eq. (jPip is verified. 

The first use we will make of the reciprocal frame is for 
the definition of two derivative operators. In fiat space 
we define the vector derivative 

V = (t"9„. (D6) 

It will be convenient, sometimes, to use vector derivatives 
in subspaces of 5D space; these will be denoted by an 
upper index before the V and the particular index used 
determines the subspace to which the derivative applies; 
For instance — a"^dm = cr^di + cr^d2 + a^d^. In 
5-dimensional space it will be useful to split the vector 
derivative into its time and 4-dimensional parts 

V = -aodt + (T'd, = -aadt + V. (D7) 

Consistently with the boldface notation for 3-dimensional 
vectors "V will be denoted by V. We will use over ar- 
rows, when necessary, to imply that the vector derivative 
is applied to a function which is not immediately on its 
right; for instance in B and in the derivative 

operator is applied to function B. 

The second derivative operator is called covariant 
derivative, sometimes also designated by Dirac operator, 
and it is defined with recourse to the reciprocal frame g" 

D = g"9a. (D8) 

Taking into account the definition of the reciprocal frame 
(jDl|l we see that the covariant derivative is also a vec- 
tor. In cases where there is a refractive index, it will be 
possible to define both derivatives in the same space. 



Vector derivatives can also be left or right multiplied 
with other vectors or multivectors. For instance, when V 
is multiplied by vector a on the right the result comprises 
scalar and bivector terms Va = V • a -I- VAa. The scalar 
part can be immediately associated with the divergence 
and the bivector part is called the exterior derivative; in 
the particular case of Euclidean 3-dimensional space it 
is possible to define the curl of a vector by curl(a) — 
-cri23VAa. 

We define also second order differential operators, des- 
ignated Laplacian and covariant Laplacian respectively, 
resulting from the inner product of one derivative opera- 
tor by itself. The square of a vector is always a scalar and 
the vector derivative is no exception, so the Laplacian is 
a scalar operator, which consequently acts separately in 
each component of a multivector. For 4 + 1 flat space it 
is 

One sees immediately that a 4-dimensional wave equation 
is obtained zeroing the Laplacian of some function 

V^* = (^-^ + V^j * = 0. (DIO) 

This procedure was used in Ref. [4] for the derivation 
of special relativity and extended in Ref. 6] to general 
curved spaces. 

APPENDIX E: COMMUTATION RELATIONS 

We examine here the commutation of operators with 
V. First of all we let us expand xAV 

xAV = {x^d2''X^di)a^^ + {x^d3'^x^d2)cr^^ + 

+ {x^di-x^d3)a^\ (El) 

The angular momentum operator associated with the 
plane of 3D bivector B is defined as 

iB=iB-(xAV). (E2) 

For instance, if we are interested in the angular momen- 
tum relative to the a'^ direction 

L^i2 = iCTi2.(xAV) = iix'^di - x^d2). (E3) 

The angular momentum operator does not commute with 
V; following Ref. [l| we have 

[B-(xAV),V] = -^B-(^AV) 

= B X V. (E4) 

Since B x V = (BV — VB)/2 we can define an operator 
B-(xAV) — B/2 which commutes with both V and T). 
The conserved total angular momentum operator is then 

Jb=Lu-^ iB. (E5) 
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We simplify the notation for the case of bivectors nor- For the K, operator we have then 
mal to frame vectors by writing = J„no^ with ct™ = 

If we ignore second derivatives (xA V)V is zero but [^^"(1 - xA V), V] = 2cr4"V - cr^° V(xA V) = 0. (E7) 

V(xAV) = 2V. (E6) 
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